www.gusucode.com > 基于matlab编程蒙特卡罗仿真设计源码程序 > 基于matlab编程蒙特卡罗仿真设计源码程序/第五讲 蒙特卡洛方法的应用/ybjzf.m
function result=ybjzf(a,b,m,mm) %a是积分的下限 %b是积分的上限 %积分函数cos(x)+2 %mm 是随机实验次数 sum=0; xrandnum = unifrnd(a,b,1,mm); for ii=1:mm sum=sum+cos(xrandnum(1,ii))+2; end result=sum*(b-a)/mm